Classical Simulation of Squeezed Vacuum in Optical Waveguide Arrays 
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We reveal that classical light diffraction in arrays of specially modulated coupled optical waveg- 
uides can simulate the quantum process of two-mode squeezing in nonlinear media, with the waveg- 
uide mode amplitudes corresponding the signal and idler photon numbers. The whole Fock space is 
mapped by a set of arrays, where each array represents the states with a fixed difference between 
the signal and idler photon numbers. We demonstrate a critical transition from photon number 
growth to Bloch oscillations with periodical revivals of an arbitrary input state, associated with an 
increase of the effective phase mismatch between the pump and the squeezed photons. 
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I. INTRODUCTION 



II. TWO-MODE SQUEEZING IN NONLINEAR 
MEDIA WITH ARBITRARY PHASE MISMATCH 
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Arrays of coupled optical waveguides enable flexible 
manipulation of optical beams [3, HI- Remarkably, the 
beam evolution can emulate the fundamental wave phe- 
nomena in different physical systems y , including the 
famous examples of Bloch oscillations ji -Q , Zener tun- 
neling , and dynamic localization [sj-flC | originally pre- 
dicted for electrons in crystalline potentials. Recently, 
it was shown that classical beam propagation in linear 
arrays can simulate quantum photonics phenomena in- 
cluding quantum walks 11| , transformation of quantum 
coherent and displaced Fock states [H|-[l5j], and gener- 
ation of photon pair s through spontaneous parametric 
down-conversion [16j . These analogies on the one hand 
enable direct visualization of the simulated photon states, 
and on the other hand suggest novel concepts for the de- 
sign of waveguide arrays offering new opportunities for 
the manipulation of optical beams. 

In this paper, we develop a special design of waveg- 
uide arrays, in which linear light propagation can sim- 
ulate the dynamics of squeezing. Optical squeezing is 
a fundamentally important quantum process, which can 
occur in media with quadratic nonlinearity. Squeezed 
light can enable high-resolution measurements, beating 
the standard quantum limit [l7| . 

The paper is organized as follows. In Sec. UH we 
overview the theory of squeezing and present the solution 
in an operator form. In the following Sec. IIII1 we for- 
mulate equations for photon numbers in Fock represen- 
tation, and demonstrate their mapping to classical light 
propagation in specially designed arrays of optical waveg- 
uides. In Sec. IIV1 we discuss the photon state dynamics 
and the critical transition to Bloch oscillations associated 
with the variation of the phase mismatch. In Sec. |V] we 
analyze the effect of periodic reversal of the sign of non- 
linear coefficient, leading to quasi-phase-matching and 
quasi-Bloch oscillations. Finally, we present conclusions 
and outlook in Sec. IVII 



We consider two-mode squeezing in the traveling-wave 
configuration. The theory of such process is well devel- 
oped [H,[l9|. In this section, we use the established ap- 
proach to derive the photon state evolution in the general 
case of varying phase mismatch between the pump and 
the squeezed photons. Such analysis is essential for our 
study, as we later demonstrate the critical transition of 
light dynamics associated with the variation of the phase 
mismatch. Specifically, we consider the Hamiltonian in 
the undepleted pump approximation, 

H = I3 s (z)a\a s + ^(z)ataj + ^j{z)a\a\ + 7*(z)a s a i , (1) 

where z is the propagation direction, a\ , t and a s ^ are the 
creation and annihilation operators for the signal and 
idler photons, 7 is proportional to the pump amplitude 
and quadratic nonlinear susceptibility, and /3 Si i charac- 
terize the phase mismatch between the signal, idler, and 
the pump. Importantly, we take into account arbitrary 
dependance of the mismatch and nonlinearity coefficients 
on the propagation coordinate. Such dependence can be 
realized experimentally, for example through domain re- 
versal [201 ] and waveguide tapering. 



A. General solution for the operator evolution 

Following [3 [l9| we obtain Heisenberg equations for 
the evolution of the operators, 



da. 



a,,H 



dz 



a,,H 



(2) 



Substituting Eq. (TTJ) into Eq. @, we obtain 

i^l = f3 s (z)a s + j(z)aj, 
dz 

dz 



(3) 
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We seek solution of Eq. (|3|) in the form 

a s (z) = U(z)a s (0) + V(z)al(0) exp[-iS(z)], 
a,{z) = [U(z)ai(0) + V(«)aj(0)] exp\iS(z)\, 

with the initial conditions 

U(0) = 1, V(0) = 0. 

After substituting Eq. (j4]) into Eq. we get 

i^E=P( z )U + 1 (z)V*, 
clz 

.dV 



and 

where 

P{z) 



= P{z)V + y(z)U* , 



5= I 6(0%, 

'0 



P,{z) + fr(z) = p 3 (z)-Pj(z) 



(4) 

(5) 

(6) 
(7) 

(8) 
(9) 



It can be checked that (|C/| 2 — | V" | 2 ) is a conserved quan- 
tity of Eqs. (HI),©. Then, we can use the representation 

U = U a exp(i(p), V — V a exp(i(p), (10) 

where 

U a = cosh(|a|), V a = — exp [i arg(a)] smh(|a|), (11) 
and 

, -i ( y \ \ ( y\~\ 

a = tanh — exp % arg 

V U J P [ H U)\ ' (12) 

ip = arg (U) . 
The operator solution can then be written as: 



a s (z) 







-i(p—i5 



u a v a 
u* V* 



a,(o) 

4(0) 
(13) 



where the two-mode squeezing or two-photon displace- 
ment operator is 



S a — exp(—iH a ) — exp(a* a s a,i — aa\a\) 



(15) 



Now considering the case of arbitrary detuning and 
nonlinearity parameters, and taking into account the 
equivalent representations discussed above, we can write 
the operator solution (TT^]) as 



a s (z) 



; exp 



K')-(^"*)«.(°?!o! 

S a exp (-iff w ) exp (-iff«) 



where 



ffW = (6 - y)a\a s , H® = (-6 - <p)o\t 



(16) 



(17) 



Accordingly, in an interaction picture solution can be rep- 
resented as 

|*(z)) = exp exp (-iff^) £ Q |*(0)) (18) 

We can now calculate the average photon num- 
ber evolution. Note that the terms exp(— iff'*)) and 
cxpHff( s )) can be neglected in this calculation, since 
these operators preserve the photon numbers: 

(n s (z)) = (al(z)a s (z)) 

= (*(o)|5£a!(o)a.(o)&|*(o)) 

= (n s (0)) + |^| 2 [l + (n s (0)) + (n l (0))] 

+ (#(0)| [u*Val(0)al(0) + UV*a s {0)ai(0)\ |*(0)) 

(19) 

If the input state has a fixed photon number, i.e. 
|tf(0)> = |ni 0) ,n[ 0) ), then 



(n,(z)}=n<°> + \V\ 2 [l+n^+n 



,<°) 



{n l {z))=nf ) + \V\ 



1 + ni n) + n 



(20) 



B. Solution in terms of standard two-mode 
squeezing operator 

Let us for a moment consider the case of zero detuning 
(Ps = Pi = 0) and constant nonlinearity (7 = const). 
Then tp = 8 = p = and the operator solution is given 
by Eq. (JI3J) with {U, V} = {U, V} a and a = i-yz. On 
the other hand, the operator solution can be equivalently 
written in a general form as (l7| 



a s (z) 



fit f a *(°) 

' " 1 4(0) 



(14) 



III. WAVEGUIDE ARRAYS DESIGN FOR 
SIMULATING FOCK STATE EVOLUTION 

A. Coupled-mode equations for the photon 
numbers 

Two-mode squeezed state in the number (Fock) basis 
can be written as 1171: 



+00 +00 



(21) 



n B =0 rii=0 
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Equations for the evolution of the state vector are 



dz 



H^(z). 



(22) 



Since the photons are generated in pairs, only the 
states with 



N = const 



(23) 



will be coupled. Since the signal and idler equations are 
symmetric, with no loss of generality we consider N > 0, 
and obtain the coupled-mode equations with modulated 
coupling for ip n = if) n<n+N , 



i—— = c-yipi + N/i^o, 
dz 

i—j— = C n 1p n -i + C n+1 lf) n+ i 

+ [n& + (n + N)pi] ip n , 



(24) 



(25) 



where n > 1 and c„ are the modulated coupling coeffi- 
cients: 



c n (z) = \Jn{n + N)j(z), 



(26) 



We now simplify the coupled equations. First, we make 
a transformation of the propagation coordinate, 



N=0 



N=l 



% ^^^^^ 



o 



2 3 



FIG. 1. (Color online) Schematics of waveguide arrays with 
varying waveguide separation to achieve the coupling \C n \ ac- 
cording to Eq. (1301) . The coupling and array design depends 
on the difference between the signal and idler photon num- 
bers: (a) N = 0, (b) N = 1, (c) N = 2. 



Z(z)=% 1 / l7(0K, 
Jo 

where 70 is a positive scaling coefficient, and obtain 
^ = ^1 + ^(70/17(^)1)^0, 

dZ 



(27) 



(28) 



[n(3 s + (n + N)(3 t ] ( 7o /| 7 (Z)|)^, (29) 



where 



where 



dZj \%Z)[ 



(35) 



To perform an optical simulation of Fock-space dynam- 
ics, one could use straight arrays of coupled waveguides, 
with parameters designed to match Eqs. ([33|) . (|34|) . This 
would be similar to Ref. [15j . but with different waveg- 
uide coupling. 



C n {z) = v/n(n + 7V) 7o exp [i q(Z)] , %Z) = 7 (z), 



and 

q(Z) = arg [j(Z)} 
Next, we make a substitution 



(30) 
(31) 



MZ) = E n (Z) exp [iqn - iN jf J|_deJ , 

(32) 

and obtain equations with positive Z-independent cou- 
pling coefficients 



.dE„ 



(33) 



dZ 



|C n |£ n _i + \C n+1 \E n+1 + p{Z)nE n , (34) 



IV. BLOCH OSCILLATIONS OF SQUEEZED 
LIGHT 



A. Transition from non-oscillatory regime to Bloch 
oscillations 



We now consider the constant (z- independent) coeffi- 
cients /3 Sj j and 7 . Then, the solution of Eqs. ©,(0) for 
the operator amplitudes can be written explicitly as 



where 



i/3 

U = cosh(6z) sinh(&z), 



7 

V = —i— sinh(foz). 



(36) 



(37) 
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We then use Eq. (|20|) to find the average photon number 
evolution for an input state with a fixed photon number, 

l*(0)> 



(n s {z)) = (ni(z)) 



= n<°> 



N 

2 



j^L|sinh(6z)| 2 



1 + N + 2n ( s 0) 



(38) 



,(0) 



where according to Eq. (f23|) . N = n\ ' — ni 

We see that the type of evolution fundamentally 

changes when the detuning crosses the critical value /3 cr . 
For detunings smaller than the critical value, |/3| < 

/3 cr , b is real, and at large z the values of U and V will 

rapidly grow corresponding to quickly increasing photon 

numbers according to Eq. (|19|) . 

For detunings larger than the critical value, \j3\ > /3 CT , 

b is imaginary, and solution becomes oscillating with the 

period 



Az 



\b\ JW^M 



(39) 




FIG. 2. (Color online) Evolution of light intensity in waveg- 
uide arrays simulating the photon number distribution in 
squeezed vacuum states, for different values of the phase mis- 
match ft: (a) 0, (b) 1, (c) 1.3, and (d) 1.86. Dashed lines show 
the average signal photon number. For all the plots, 7 = 1 
and N = 0. 



At distances z — mAz, where m — 1, 2, 3, 4, . . ., the so- 
lution returns to the input state. Note that such dy- 
namics will occur for any input, in lattices corresponding 
to any value of N. Considering the realization in pho- 
tonic lattices, we can call this regime Block oscillations 
of squeezed light as a generalization of spatial Bloch os- 
cillations in waveguide arrays 04^1 or squeezed quantum 
bouncing ball as a generalization of quantum bouncing 
ball in photonic lattices [U H2| • 

In the special case of equal signal and idler photon 
numbers, when N — 0, the coupling and detuning coeffi- 
cients in Eqs. (|3"3")l . (fM]) directly correspond to the previ- 
ously studied lattice with linearly increasing hopping rate 
and on-site potential (23|. In Ref. [23j . a transition be- 
tween Bloch oscillations and non-periodic dynamics was 
first identified, which agrees with our findings for a more 
general coupling dependence corresponding to arbitrary 
N. Actually, we find that the lattice modulation param- 
eters in Eqs. (|3"3")l . (13~4"|) belong to a previously identified 
class of lattices with commensurate energy levels [24| . 
where Bloch oscillations were predicted. Specifically, our 
model corresponds, up to a gauge transformation, to 

(1),(2) in Ref. [II with Fi(ra) = ^ /2 q(n + N + 1) 



Eqs. 



and Fzin) = 7 q n, where q is defined to satisfy the 
relation (q + q^ 1 ) = —pj 1 ^ 2 - 

There were also recent studies of Bloch-like oscillations 
in Jaynes-Cummings [2f| and Glauber- Fock photonic lat- 
tices [3, [H| , where the oscillation period was found to be 
inversely proportional to the linear detuning between the 
waveguides (/3 in our notations). These results can be ob- 
tained as a limiting case of our study by taking N — > oo 
and 7 = 70 — >• joo/\/N- In this case, C n — > "fooVn and 
Az = 2tv/\/3\, as for the Glauber-Fock lattice considered 
previously [hi. Il5|. 



B. Evolution of the input vacuum state 

We first illustrate the general results for the squeezed 
vacuum state. The vacuum state at the input corre- 
sponds to the initial conditions 

N = 0, ^(z = 0) = l, ^ n>1 (z = 0) = 0. (40) 

Combining the expression for the squeezed vacuum 
state [13 with Eq. CDS), we obtain 



2itpn 



cosh(|a|) 



{-l) n e me [tanh(|a|)]™ (41) 



For convenience, we can choose the scaling coefficient 
7o = ItIj such that Z = z, and then the normalized 
mode amplitudes are found as 



E n (Z) = [cosh 2 (bZ) 



■ P 2 \b~ l siDh(bZ)\ 2 ] 



-1/2 



-«|7|6- 1 sinh(6Z) 



cosh(6Z) + i$b~ x sinh(6Z) 



(42) 



We plot the characteristic dependencies of the mode 
intensities in Figs.[5Ja-d) for different values of the phase 
mismatch. In the phase-matched regime, the photon 
number increases rapidly, as marked with the dashed line 
in Fig.[2ja). For phase mismatch chosen at the threshold, 
P — f3 C r, the photon number also grows continuously, but 
at a slower rate, see Fig.^b). However, as the mismatch 
is increased above the threshold value, the photon num- 
ber growth becomes bounded and the photon distribution 
periodically returns to the input state, see Figs. [UJc,d). 
In particular, the oscillation period is Az — 2 for the 
mismatch parameter in Fig. [U[d). 
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FIG. 3. (Color online) Synchronous Bloch oscillations in 
waveguide arrays corresponding to various input conditions: 
(a,c,e) input at the waveguide n = corresponding to — 
and (b,d,f) input at the waveguide corresponding to ni = 1. 
Dashed lines show the average signal photon number. The 
photon population differences between the signal and idler 
photons are (a,b) N = 0, (c,d) N = 1, (c,d) N = 2, and 
these are simulated using arrays with different couplings as 
illustrated in Fig. [T] For all the plots, 7=1 and (3 = 1.86. 



C. Synchronous Bloch oscillations of squeezed 
quantum bouncing ball 

We now illustrate that for arbitrary input, the Bloch 
oscillations occur synchronously, with the same oscilla- 
tion period according to Eq. (|39|) . We fix the values of the 
nonlinear and phase mismatch coefficients, and perform 
numerical simulations corresponding to different initial 
conditions in Fig. [3J The three rows in this figure show 
the evolution for the photon number differences between 
the signal and idler modes N = 0,1,2, and these corre- 
spond to accordingly different configurations of waveg- 
uide arrays as illustrated in Fig. [TJ The left and right 
columns in Fig. [3] correspond to different excitation con- 
dition of each array, where the light is coupled to the first 
or second waveguide at the input, respectively. Such cou- 
pling corresponds to the input states with zero (n^ = 0) 
or one = 1) photon in the signal mode at the in- 

put. We observe that the photon number distributions 
are strongly dependent on the input number of photons 
in the signal and idler modes. Nevertheless, in agreement 
with the analytical predictions, the full revival of an in- 
put state is observed after on Bloch oscillation period. 

Such features of Bloch oscillations with the periodic 



revivals indicate that the spectrum of eigenmodes in the 
lattices is equidistant, forming a semi- infinite Wannicr- 
Stark ladder with the spacing between the levels of 
2ir/ Az. Moreover, this spectrum is the same for different 
arrays corresponding to various values of N. Therefore, 
we have identified an infinite family of lattices with dif- 
ferent non-uniform coupling parameters, but the same 
equidistant spectra. This can offer new opportunities for 
optical beam manipulation. 



V. PERIODIC RECONSTRUCTION AND 
SQUEEZING IN QUASI-PHASE-MATCHED 
STRUCTURES 



In media with quadratic nonlinearity, the efficiency 
of interactions can be increased under the presence 
of phase mismatch through the quasi-phase-matching 
(QPM) technique [20| . With this approach, the sign 
of the quadratic nonlinear coefficient is reversed along 
the propagation direction. We analyze the effect of such 
modulation of the simulation of squeezing dynamics in 
waveguide arrays, considering a periodic modulation of 
the sign of the nonlinear coefficient. Specifically, we take 
7(as) = 70 > for < x < L/2 and 7(35) = -7(0; - L/2) 
for x > L/2, where L is the modulation period. In terms 
of the coupled-mode Eqs. (|3"3"]) . ([55)1 . each the sign rever- 
sals correspond to p(Z) = K^ m 5(Z — mi/2), where 
S(Z) is Dirac delta function. Accordingly, at these posi- 
tions the phase of the wave function is abruptly changed 
as E n (mL/2 + £) = (-l) n E n (mL/2 - £) for £ -> +0. 
Such phase modulation was suggest ed p reviously for the 
realization of image reconstruction 26, 27], however we 
find that in our system the possibility of input state re- 
construction strongly depends on the phase mismatch. 

We present the simulation results in Figs. QJa-d) for 
the evolution of an input vacuum state with different 
values of the phase mismatch. In the case of exact phase 
matching (/J = 0), the reversal of the sign of nonlinear 
coefficient effectively reverses the propagation dynamics, 
such that after each modulation period the input state is 
restored, similar to Rcfs. [2(|[27|. As a phase mismatch is 
introduced, the reversals are suppressed, see Figs.0Jb,c). 
This agrees with the principles of QPM, that the effective 
phase mismatch is modified by the reciprocal wavevector 
2tt I ' L. The fastest increase of the photon number is visi- 
ble in Fig.[3Jc), when the mismatch (f3 — 1.86) is chosen 
such that the Bloch oscillation period in the absence of 
QPM [Fig. [2(d)] exactly matches the QPM period. For 
particular values of mismatches, the revival of the input 
state can occur as shown in Fig. 0Jd). The latter case 
resembles the regime of quasi-Bloch oscillations and dy- 
namic localization 1281 . 
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FIG. 4. (Color online) Effect of periodic sign reversal of the 
nonlinear coefficient with the period L — 2 for different values 
of the phase mismatch ft: (a) 0, (b) 1, (c) 1.86, and (d) 2.5. 
Dashed lines show the average signal photon number. For all 
the plots, |7(x)| = 1 and N — 0. 



directly model the photon number evolution in the pro- 
cess of quantum two-mode squeezing, in the general case 
of coordinate-dependant phase mismatch and nonlinear- 
ity coefficients. We identified a phase-mismatch depen- 
dent transition between a gradual beam diffraction corre- 
sponding to an increase of the squeezed photon numbers, 
and Bloch oscillations with the periodic revival of an ar- 
bitrary input state. Additionally, we have illustrated a 
different regime of periodic revivals associated with the 
periodic reversal of the effective nonlinear coefficient, re- 
sembling the regimes of input reconstruction and quasi- 
Bloch oscillations. 



We anticipate that our results on the families of waveg- 
uide arrays with modulated parameters may stimulate 
further studies towards novel opportunities of employing 
such lattices for manipulation of optical beams including 
all-optical switching of intense beams @, |j| and control of 
spatial entanglement through quantum walks with corre- 
lated photons [HI. 



VI. CONCLUSIONS 

In conclusion, we have demonstrated that arrays of op- 
tical waveguides with specially modulated coupling can 
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